1 



Anomalous diffusion of a particle in an aging medium 

Noelle POTTIER 

Federation de Recherche CNRS 2438 "Matiere et Systemes Complexes" 

and 

Groupe de Physique des Sohdes, CNRS UMR 7588, Universites Paris 6 et Paris 7, 
2, place Jussieu, 75251 Paris Cedex 05, France 

and 

Alain MAUGER 

Laboratoire des Milieux Desordonnes et Heterogenes, CNRS UMR 7603, Universite Paris 6, 
Tour 22, Case 86, 4 place Jussieu, 75252 Paris Cedex 05, France. 



Abstract 

We report new results about the anomalous diffusion of a particle in an aging 
medium. For each given age, the quasi-stationary particle velocity is governed by 
a generalized Langevin equation with a frequency-dependent friction coefficient 
proportional to |ci;|'^~^ at small frequencies, with < 5 < 2. The aging properties 
of the medium are encoded in a frequency dependent effective temperature T^s, (iv). 
The latter is modelized by a function proportional to |c<;|° at small frequencies, with 
a < 0, thus allowing for the medium to have a density of slow modes proportionally 
larger than in a thermal bath. Using asymptotic Fourier analysis, we obtain the 
behaviour at large times of the velocity correlation function and of the mean square 
displacement. As a result, the anomalous diffusion exponent in the aging medium 
appears to be linked, not only to S as it would be the case in a thermal bath, but 
also to the exponent a characterizing the density of slow modes. 
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The question we want to address in the present paper is that of the diffusion 
(possibly anomalous) of a particle in an out of equilibrium environment. Generally 
speaking, the medium in which a diflFusing particle evolves may be, or not, in a state 
of thermodynamic equilibrium. For instance, when it is constituted by an aging 
medium such as a glassy colloidal suspension of Laponite [l]-[3], the environment 
of a diffusing particle is not in thermal equilibrium. This renders the analysis of the 
particle motion much more involved, as compared with the thermal equilibrium 
case. 

As well-known, the motion of a diffusing particle in a stationary medium is 
properly described by a generalized Langevin equation, in which the effect of the 
environment is encoded in a friction term and a noise term. When the environment 
is a thermal bath, the noise spectral density and the real part of the friction 
coefficient are related by the fluctuation-dissipation theorem of the second kind 
[4], [5]. Thermalizing particle variables such as the particle velocity then satisfy the 
fluctuation-dissipation theorem of the first kind i.e. the Einstein relation [4], [5]. 
In both theorems the temperature of the bath plays the central role. 

None of these theorems is valid when the environment ages, since then it 
is out of equilibrium. It has been proposed to associate to an aging system a 
violation factor of the fluctuation-dissipation theorem, or an effective temperature 
[6], [7]. This latter quantity is considered as a function TesX^jiw) of the angular 
frequency u of an external perturbation and of the age of the system. 

In the present paper we analyze in details how this concept allows for the 
description of diffusion in an aging medium. We propose to modelize TQff_{u>,t^) 
by an inverse power-law of lv at small uj, in accordance with the fact that, in an 
aging medium, the density of slow modes is proportionally enhanced as compared 
to what it is in a thermal bath (for the same real part of the friction coefficient). 
We show that this modelization leads to an anomalous diffusion exponent which 
depends, not only on the exponent characterizing the real part of the friction 
coefficient, as it is the case at equilibrium, but also on the exponent characterizing 
Teff.. This constitutes the main new result of the paper. 

The paper is organized as follows. In Section 2, we briefly recall the main 
features of the diffusion of a particle in a stationary medium. In particular, we 
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carefully separate out the properties valid at thermal equilibrium (namely, the 
two fluctuation-dissipation theorems and the regression theorem) from those which 
only require stationarity. In Section 3, we consider the general situation of a par- 
ticle diffusing in an out of equilibrium environment. Assuming a bilinear coupling 
between the particle and the environment, we derive an equation of motion for the 
particle. This equation involves an environment response function and a fluctuat- 
ing force. In contrast to the equilibrium case, there is no speciflc relation between 
these two quantities. In Section 4, we show on general grounds how, when there 
exists two well separated time scales respectively characterizing the times perti- 
nent for the measuring processes and the waiting time or the age of the system, a 
quasi-stationary regime can be deflned, as pictured by frequency-dependent gen- 
eralized susceptibilities (parametrized by tu,). We then focus the study on the 
environment response function and on the particle velocity response function. We 
show how, in the quasi-stationary regime which takes place when ujtyj ^ 1, an 
effective temperature Teff.(a;, t^j) can be deflned, which modifles in a consistent 
way both the Einstein relation and the relation between the noise spectral density 
and the real part of the friction coefficient. In Section 5, we propose a modeliza- 
tion of Teff. by an inverse power-law of uj at small a;, which corresponds to an 
enhanced density of slow modes with respect to equilibrium. Using asymptotic 
Fourier analysis, we compute the anomalous diffusion exponent. As a result, it 
depends on both exponents characterizing the real part of the friction coefficient 
and the effective temperature. This result shows that independent measurements 
of the real part of the particle mobility and of the mean square displacement will 
in turn allow for the determination of the exponent associated with the small oj 
behaviour of Teff ., which describes the density of slow modes in the aging medium. 
Finally, Section 6 contains our conclusions. 



2. Difrusion in a stationary medium 

To begin with, we recall in this Section some known results about the diffusion 
(possibly anomalous) of a particle in a stationary medium. For further purpose, 
we will separate out the results valid for any stationary medium from those which 
are only valid for a thermal bath. 

The motion of a diffusing particle of mass m evolving in a stationary medium 



4 



is usually described by the generalized Langevin equation [4] , [5] , 

m^ = -mj 7(t - t') v{t') dt' + F{t), ^ = (2-1) 

in which F{t) is the Langevin random force acting on the particle and ^{t) is a 
retarded friction kernel. In Eq. (2.1) it is assumed that the diffusing particle and 
the surrounding medium have been put in contact in an infinitely remote past, as 
pictured by the lower integration bound — oo in the retarded friction term. Both 
F{t) and the solution v{t) of the generalized Langevin equation (2.1) can be viewed 
as stationary random processes. Their spectral densities are linked by 

C^,{u;)^\li{u;)fCFF{u;), (2.2) 

where 

fiiu) = — (2.3) 

denotes the frequency-dependent particle mobility (7(0;) is the Fourier transform 
of 7(t), as defined by 7(0;) = ^{t) e'^^ dt). 

2.1. Relaxation of the average velocity 

As shown by the expression (2.3) of the mobility, the function 'y{uj) character- 
izes the relaxation of the average particle velocity. For instance, with a power-law 
behaviour of 7(0;) (namely 3?e7(c<;) oc with < 6 < 2), one has [8]: 

{v{t)) = v{t = 0) E2-5 [-M^-'] . (2.4) 

In Eq. (2.4), Ea{x) denotes the Mittag-Leffler function -"^ of index a [9], [10], and 
the 5-dependent frequency cos acts as an inverse relaxation time. For S = 1, 



^ The Mittag-Leffler function is defined by the series expansion 

Ea{x) = V — — --, a > 0, 

^ T{an + 1) 

n=0 

where F is the Euler Gamma function. The Mittag-Leffler function Ea{x) reduces 
to the exponential when a = 1. The asymptotic behaviour at large x of the 
Mittag-Leffler function Ea{x) is as follows: 

E^{x) - ~ ^ x>l. 
X L [1 — a) 
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the function £^2-<5 [— ('^5^)^ ^] reduces to a decreasing exponential {Ei{—uit) = 
e"'^^*), while, for all other values of 5, it decays algebraically at large times 

{E2-s[-{costy-^] ^ {u;sty-yr{S-l)) [10]. 

2.2. Velocity correlation function and mean square displacement 

Applying the Wiener-Khintchine theorem, one obtains the velocity correlation 
function as the inverse Fourier transform of Cyy{Lo), that is, in terms of the noise 
spectral density CpFi'-^)- 

{vitHt')) r ^ e-^-(*-*') CMo^) ^J— . ■ (2.5) 

Eq. (2.5) can be rewritten as: 

U(t)v(t')) ^ 1- r '^^-Mt-t')CFF{oji r 1 1 1 

w? J_^2tt 2 3?e7(a;) \7(a;)-za; 7*(a;)+ia;J ■ 

(2.6) 

From the velocity correlation function, one can compute the mean square displace- 
ment of the diffusing particle (a quantity which cannot be written under the form 
of a Fourier integral) as a double integral over time: 



/\x^{t) = {[x{t) - x{0)f) = 2 [ dU r dt2 {v{ti)v{t2)). (2.7) 

Jo Jo 

2.3. Case of a thermal bath: fluctuation-dissipation theorems 

Let us now consider the particular case of diffusion in a thermal bath. The 
linear response theory, applied to the particle velocity considered as a dynamical 
variable of the isolated particle-plus-bath system, allows to express the mobility 
in terms of the equilibrium velocity correlation function. Since the mobility /i(a;) 
is nothing but the generalized susceptibility Xvxi^^), one has 

^^"^^^kf i (^(^M^^"'^^' (2-8) 

where T is the bath temperature [4], [5]. The result (2.8) constitutes the first 
fiuctuation-dissipation theorem (FDT). Accordingly, the velocity spectral density 
is given by: 

/oo 
{v(t)v)e"^^ dt = kT21ken{uj). (2.9) 
-oo 
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Introducing the frequency-dependent diffusion coefficient D{u) as defined by 

nOO 

D{uj)= {v{t)v) cos ojtdt, (2.10) 
Jo 

one can rewrite Eq. (2.9) as the well-known Einstein relation between the difi^usion 
coefficient and the real part of the mobility: 

^("^ =kT. (2.11) 



3?e /x(a;) 



Since //(a;) can be extended into an analytic function in the upper complex half 
plane {Qmoj > 0), Eq. (2.8) on the one hand, and Eqs. (2.9) or (2.11) on the 
other hand, constitue two equivalent formulations of the first FDT. 

Then, using the expression of C„^,(a;) as given by the first FDT (2.9), that is 

C„.(.) = ^^^!?pH , (2.12) 

one gets from Eq. (2.2) the noise spectral density, 

Cff{co) = mkT2^e^{u;), (2.13) 
a relation from which one deduces the expression of 7(0;): 

^^""^^mi^ i (^(t)^)'^'^'dt. (2.14) 

Eq. (2.14) is known as the second FDT [4], [5]. Eq. (2.13) on the one hand, and 
Eq. (2.14) on the other hand, are two equivalent formulations of the second FDT. 
The second FDT can also be established directly by applying the linear response 
theory to the force exerted by the bath on the particle, this force being considered 
as a dynamical variable of the isolated particle-plus-bath system. We will come 
back to this point in Section 3. 

Thus, when the particle environment is a thermal bath, the two fluctuation- 
dissipation theorems are valid. In these two theorems the bath temperature T 
plays an essential role. Under its form (2.9) or (2.11), the first FDT involves the 
spectral density of a dynamical variable linked to the particle (namely its velocity), 
while, under its form (2.13), the second FDT involves the spectral density of the 
random force, which is a dynamical variable of the bath. 
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2.4. Regression theorem 

When the particle environment is a thermal bath, we may introduce in the 
expression (2.6) of {v{t)v{t')) the value of the ratio CFF{i^)/'^^e^{uj) as given by 
the second FDT (2.13). One gets: 

{vit)v{t')) ^ — r |^e--(*-*') i , ! . + , ] . ]. (2.15) 

In the above integral, the function l/[^{u!) — iuj] is analytic for ^mu > 0, while 
the function + iuj] is analytic for Qmuj < [5]. Computing {v{t)v{t')) 

ioT t — t' > 0, one chooses as an integration contour the semi-circle of large radius 
in the lower complex half plane. The second term does not contribute. One thus 
has: 

{v{t)v{t')) = ^ r ^ t-t'>0. (2.16) 

Formula (2.16) displays the fact that, when diffusion takes place in a thermal bath, 
the velocity correlation function is characterized by the same law as the average 
velocity [8]. This result constitutes the regression theorem, valid at equilibrium 
for any 7(0;). 

In particular, with the model for 7(0;) considered before (3^67(0;) oc |a;|'^~-'^ 
with < 5 < 2), one has: 

kT 

{v{t)v{t')) = —E2-5[-Mt-tf-']. (2.17) 
In the particular case t = t' , one gets from Eq. (2.17) the equipartition result : 

(.^> = ^. (2.18) 

3. Out of equilibrium environment 

The fully general situation of a particle diffusing in an out of equilibrium 
environment is much more difficult to describe. Except for the particular case of 
a stationary environment, the motion of the diffusing particle cannot be described 
by the generalized Langevin equation. A more general equation has to be used. 
The fluctuation-dissipation theorems are a fortiori not valid. However, one can try 
to extend these relations with the help of an age and frequency dependent effective 
temperature, such as proposed and discussed for instance in [6], [7]. 
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An equation of motion will be derived in the present Section for the case of 
a bilinear coupling between the particle and the environment. As for the out of 
equilibrium extension of the linear response theory, it will be discussed in the next 
Section. 

3.1. Equation of motion of a particle linearly coupled to an out of 
equilibrium environment 

The validity of the generalized Langcvin equation (2.1) is restricted to the case 
of a stationary medium. In other cases, for instance when the diffusing particle 
evolves in an aging medium such as a glassy colloidal suspension of Laponite [1]- 
[3], another equation of motion has to be used. We shall derive such an equation 
by assuming a bilinear coupling, of the form — $x, between the particle and its 
environment. With such a coupling, the particle equation of motion reads 

m^ = m, (3.1) 

where $(t) denotes the global force exerted by the environment. In turn, the 
latter is perturbed by the coupling with the particle. The linear response relation 
yielding the average force exerted by the environment is of the form 



where the causal function X4>4>{tji') {t > t') is a linear response function of the 
surrounding medium (in Eq. (3.2), it is assumed that the particle and the envi- 
ronment have been put in contact at time t = 0). The particle equation of motion 
(3.1) takes the form 



where ^E{t) is a fluctuating force of zero mean generated by the environment. 
When the latter is out of equilibrium, there is no specific relation between the 
response function x^^{t,t') and the correlation function {^E{t)^E{t'))- In the 
most general situations, neither ^Eit) nor v{t) can be viewed as stationary random 
processes. Before considering these situations, let us come back briefly on the 
particular cases, first, of a stationary medium, second, of a thermal bath, in order 
to see how the properties described in Section 2 can be retrieved. 




(3.2) 




(3.3) 
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3.2. Stationary medium case 

Let us consider again the particular case of a particle diffusing in a stationary 
medium, in order to see how the generalized Langevin equation (2.1) can be de- 
duced from the more general Eq. (3.3). When the medium is stationary, x$$(t, t') 
reduces to a function of t — t' (x$$(t,t') = X<s>'S>{'t — t'))- Introducing then the 
causal function j{t) as defined by 

xMt) = -rn^, (3.4) 
and integrating by parts, one gets from Eq. (3.3): 

m^ = -m ^{t-t')v{t')dt' + ^E{t) + m^{0)x{t)-m^{t)x{0). (3.5) 
dt Jq 

Let us now denote by ti (instead of 0) the initial time at which the particle has 
been put in contact with its environment. Eq. (3.5) then writes: 

dv 

m-^ = -m / ^{t-t') v{t') dt' + <5Eit) +m^{U)x{t-ti) -m^{t-U)x{U). (3.6) 

Eq. (3.6) can be identified with the generalized Langevin equation (2.1) by letting 
ti — > — oo (one has then ^(ti) = 0), and assuming that the term —m^{t — ti)x{ti), 
which depends on the initial position of the particle, is negligible in this limit. 
The fluctuating force $£;(t) then identifies with the Langevin force F{t). The 
fluctuating force and the particle velocity can be viewed as stationary random 
processes. When equilibrium is not realized, the fluctuation-dissipation relations 
(2.8) and (2.14) and the regression theorem (2.16) are not valid. 

3.3. Thermal bath case 

When equilibrium is realized, the stationary medium constitutes a thermal 
bath. The linear response theory, applied to the bath dynamical variable $e(^), 
yields 

X<,^{t - t') = PG{t - t')^{F{t)F{t')), P = {kTr\ (3.7) 

that is, using Eq. (3.4) and the causality property, and the fact that at equilibrium 
^E{t) is identical with the Langevin force F{t): 

= " ^'^ (^(Oi^(^')>- (3.8) 

Eq. (3.8), through Fourier transformation, yields the second FDT (Eq. (2.14)). 
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4. Out of equilibrium linear response theory 



When the environment is not stationary, response functions such as x^^{t, t') 
and Xvx{t, t') (response function of the particle velocity) depend separately on the 
two times t and t' entering into play, and not only on their difference t = t — t'. 
However, the time difference r, also called the measure time or the observation 
time, continues to play an essential role in the description. For this reason, it 
has been proposed to define time and frequency dependent response functions 
as Fourier transforms with respect to r of the corresponding two-time quantities 
[6], [7], [11]. The time t, which represents the waiting time or the age of the system, 
then plays the role of a parameter. 

Let us first briefly recall these definitions and examine in which conditions 
the age and frequency dependent response functions share the analytic properties 
of the corresponding stationary quantities. 

4.1. Age and frequency dependent response functions 

Considering the response function x(t, t') as a function of t and of the measure 
time T = t — t' (i.e. x(t, t') = Xi(t, r)), one introduces the Fourier transform of xi 
with respect to r (for a fixed t): 



Due to the causality of x{t,t'), the lower integration bound in formula (4.1) is 0. 
As for the upper integration bound, it is equal to t minus the lower bound over 
t' . If one assumes that the perturbation is applied at = 0, the upper integration 
bound is equal to t. One then writes 




(4.1) 




(4.2) 



that is: 




(4.3) 



The Fourier relation (4.1) can be inverted, which yields: 




(4.4) 
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The time t can be interpreted as the waiting time or the age tyj of the system 
under study. In other words, one has introduced [6], [7], [11]: 

Xi{u;,t^)= x{tw,t^-r)e''^^dT. (4.5) 

^0 

Note that the stationary regime is obtained, not only by assuming that x(t, t') is 
invariant by time translation {x{t,t') — Xi{t))j but also, moreover, that the age 
of the system tends towards infinity: 

X(cc;)- lim / xi{r)e"^^dT. (4.6) 

til) >oo Jq 

However, from the point of view of linear response theory, the definitions 
(4.1) or (4.5) suffer from several drawbacks. Actually, the function xi(a;,tty) as 
defined by Eq. (4.5) is not the Fourier transform of the function Xi{tw,T), but 
a partial Fourier transform computed in the restricted time interval < r < t^. 
As a consequence, it does not possess the same analyticity properties as the gen- 
eralized susceptibility x(a;) defined by Eq. (4.6). While the latter, extented to 
complex values of u, is analytic in the upper complex half plane C^muj > 0), 
the function Xii'^^tw) is analytic in the whole complex plane. As a very simple 
example, consider the exponentially decreasing response function 

x{t,t') = @{t-t')e-^^'-''\ (4.7) 

One has^: 

Xii^.t^) = [l-e-(^-^-)*-]. (4.8) 

7 — zo; 

The corresponding generalized susceptibility is obtained by taking the limit 
tu; ^ oo in Eq. (4.8): 

X(a;) = ^— (4.9) 
7 — za; 

In this example, xi'^) given by Eq. (4.9) has a pole in the lower complex half 
plane (a; = —i'y)- This pole can be traced back to an eliminable singularity in the 
r.h.s. of Eq. (4.8), so that xi('^,^w) is analytic everywhere. This is an important 



^ Despite of the time translational invariance of x{t,t'), the partial Fourier 
transform Xi(<^) t^) effectively depends on t^j, due to the fact that the perturbation 
is applied at time and not — oo. 
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drawback, since the information about the nature of the modes of the unperturbed 
system, which is contained in the poles of the generahzed susceptibihty, is not 
contained in the partial Fourier transform Xi(a;, tw) as defined by Eq. (4.5). 

4.2. Quasi-stationary regime: introduction of an effective tempera- 
ture 

However, if utw ^ 1, one can assume that xi(w, varies very slowly with 
in the range of values of oj of interest [6] , [7] . It exists then two well separated time 
scales, respectively characterizing the times pertinent for the measuring process, 
and the waiting time or the age of the system. In this case, xii^^^w) describes 
a quasi-stationary regime taking place for a given waiting time. Considered as a 
function of u, xii'^^tw) is assumed to have the analytic properties required from 
a generalized susceptibility (i.e. analyticity in a complex half plane, namely the 
upper one with our definition of the Fourier transformation). From now on, we 
will drop out the parameter for simplicity, being it understood that the analysis 
is valid for a given tw 

Let us now come back to the specific problem of the diffusion of a particle in an 
out of equilibrium environment. In a quasi-stationary regime, the particle velocity 
obeys the generalized Langevin equation (2.1) The generalized susceptibilities of 
interest are the particle mobility //(a;) = Xvxi^^) and the friction coefficient 7(0;) = 
— (l/ma;)x4>*(<^) (this latter formula being a consequence of the relation (3.4) 
between ^{t) and x$$(t)). The results of the linear response theory as applied to 
the particle velocity, namely the first FDT (2.8) and the Einstein relation (2.11), 
are not valid out of equilibrium. One can then try to extend the linear response 
theory with the help of an effective temperature [6], [7]. More precisely, we shall 
write the relation between the velocity spectral density C„„(a;) = J^^{v{t)v)e^'^^dt 
and the real part of the mobility as 

Cyy{u;) = kTe{fXu;)2^efi{uj), (4.10) 

which amounts to assume a modified Einstein relation: 

Eqs. (4.10) and (4.11) define a frequency dependent effective temperature Teff.(a;). 
Note, however, that it is not possible to rewrite the first FDT under a form similar 
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to Eq. (2.8) with T^uXlj) as defined by Eqs. (4.10) and (4.11) in place of T: 

Nevertheless, Tefj. (a;) can consistently be used in the relation between the 
noise spectral density Cff{uj) = j'^^{F{t)F))e^^'^dt and the real part of the fric- 
tion coefficient. Indeed, from Eq. (2.2), valid in any stationary regime, and 
Eq. (4.10), one deduces: 

CFF{oj)^mkT^sXuj)2^e'^{u). (4.13) 

But, as for the first FDT, it is not possible to rewrite the second FDT under a 
form similar to Eq. (2.14) with Tefj. (a;) in place of T: 

Given 3^67(0;) and Tefj.(a;), it is possible to compute by Fourier analysis the 
out of equilibrium velocity correlation function {v{t)v{t')), and to deduce from it 
the associated mean square displacement Aa:^(t). The results can be compared 
with the ones obtained at equilibrium (Teff. ('^) = T) for the same 9?e7(a;) [8]. 
Such a program will be carried out in the following Section. 

5. Anomalous diffusion in an aging medium 

In terms of the effective temperature Teff.(a;) as defined by Eq. (4.10) and 
Eq. (4.11), one has: 

v{t)v{t')) = r ^e-^-(*-*') ^^^2Ke/i(u;). (5.1) 
J -00 27r m 

Formula (5.1) displays the fact that, when diffusion takes place in an aging medium 
as pictured by a frequency-dependent effective temperature Teff.(a;), the velocity 
correlation function is not characterized by the same law as the average veloc- 
ity. In other words, the regression theorem (2.16) is not valid. Accordingly, the 
equipartition result (2.18) is not valid either. 

To characterize at best the properties of {v{t)v{t')) in an aging medium, we 
shall take for 9^67(0;) the same function than in our equilibrium study [8], namely 
a function behaving like a power-law characterized by the exponent 5 — 1: 

Ue-f{oj) = -fs (^i^) , |a;| < Uc. (5.2) 
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We limit ourselves to the case < 5 < 2, in which an initial fluctuation of the 
particle velocity relaxes towards zero at large times, as displayed by Eq. (2.4). In 
Eq. (5.2), ujc denotes a cut-off frequency typical of the environment, and cD -C 
is a reference frequency allowing for the coupling constant rjs = rwys to have the 
dimension of a viscosity for any S [12]. Since 6 < 2, the real part of the mobility 
— l/{m[y{u>) — iu>]) varies like uj^~^ for uj <^ iv: 

j£e//(a;)~ — sm — , \u\<t^u>. (5.3) 

777,75 V cD / 2 

Provided that the behaviour of Teff.(a;) at small u is known, the Fourier 
relation (5.1) allows for a simple asymptotic analysis of the large time behaviour 
of the out of equilibrium velocity correlation function. 

5.1. Modelization of the function TeffX"^) 

Let us assume the following behaviour of Teff.(a;) at small oj: 

Teff.(a;)~T^ , |a;|<a;o, a<0. (5.4) 
\a;o / 

In Eq. (5.4), both the exponent a and the characteristic frequency ujq <^ (jj may 
depend on the waiting time (a = a{tw), coq = u!o{t^)). The frequency uq is a 
decreasing function of [13], [14]. It separates low frequencies (i.e. slow modes), 
for which the modelization (5.4) applies, from high frequencies (i.e. fast modes), 
for which one has: 

Tes.{i^) ^T, > Wo- (5.5) 

The chosen modelization of the effective temperature in an aging medium insures 
that the density of slow modes in the noise is proportionally larger than in a 
thermal bath at temperature T [14]. Note that the equilibrium situation would 
correspond to ujo = 0, since then Teff.(a;) = T for any cu. It can also formally be 
retrieved by taking the limit a — > 0~ in the low frequency modelization (5.4). 

5.2. Asymptotic analysis 

Eq. (5.2), together with the low frequency expression (5.4) of Teff.(a;), shows 
that the velocity correlation function is the Fourier transform of a function varying 
like a;'*"''^"'^ at small uj (i.e. |a;| <S a;o <S u)). The Fourier integral (5.1) converges 
provided that the condition 

a>S-2 (5.6) 
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is fulfilled. According to the usual properties of the Fourier transformation, the 
correlation function {v{t)v{t')) behaves like \t — t'\~'^'^^~'^ for large values of \t — t'\. 
More precisely, one has: 

After two integrations over time (see formula (2.7)), one obtains the large time 
expression of the mean square displacement: 

Ax\t) ~ 2 - (a;ot)-" {Cot)' — — ^ ^ , \ . (5.8) 

At this stage, it is convenient to introduce the a- and 5-dependent frequency 
ujo,,5 as defined by 

-, • {5 — Oi)TT 

In terms of uJa,Sj one has: 

Out of equilibrium, the anomalous diffusion exponent, as defined by Ax^{t) ~ t^, 
is: 

1^ = 6 -a. (5.11) 

Note that the condition (5.6) insures that < 2, as required since the motion 
must not be kinematical. 

The equilibrium value z^eq. = S of the diffusion exponent is recovered, as it 
should, by making a = in Eq. (5.8). Moreover, the full equilibrium expression 
of the mean square displacement, that is [8] 



kT 1 {ustY 
'm ^ r(5 + 1) ' 



Ax\t%^. - 2 2 wr , .N ^ ^st > 1, (5.12) 



with the (5-dependent frequency cos as defined by 

1 1 



^r'=7^-^-^' (5-13) 



sm 2 
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is also retrieved from Eq. (5.7) in this limit, since 



UJS = lim UJa,5- 



(5.14) 



5.3. Discussion 

The value (5.11) of the anomalous diffusion exponent v results from the fre- 
quency dependence of both the real part of the friction coefficient and the effective 
temperature. 

The out of equilibrium noise spectral density is given by Eq. (4.13), that is, 
with the real part of the friction coefficient as given by Eq. (5.2) and the effective 
temperature as modelized by Eq. (5.4): 



As indicated above, for a given 3?e7(a;), the negative value chosen for a signifies 
that, in the out of equilibrium medium under study, the weight of slow modes in 
the noise is increased with respect to the value it would have in a thermal bath 
at temperature T. Accordingly, the anomalous diffusion exponent is enhanced : 

From the experimental side, the diffusion exponent v can be deduced from 
mean square displacement measurements, while the value of S can be obtained 
independently from mobility measurements. From both types of results, the expo- 
nent a characterizing the low frequency dependence of Teff . (u) and thus describing 
the density of slow modes in the aging medium can be deduced. 



We have studied the anomalous diffusion of a particle in an aging medium. In a 
quasi-stationary regime, the particle velocity is governed by a generalized Langevin 
equation. The medium being out of equilibrium, the fluctuation-dissipation theo- 
rems are not valid. Both the Einstein relation and the relation between the real 
part 3?e 7(0;) of the friction coefficient and the noise spectral density can be written 
in a modified way which involves an age and frequency dependent effective tem- 
perature Teff. (a;). We have first carefully made precise the conditions under which 





(5.15) 



6. Conclusion 
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such an effective temperature can be used. Then, we have studied the velocity- 
correlation function and the mean square displacement by means of asymptotic 
Fourier analysis. 

Two main parameters enter into play, d and a, which characterize the low fre- 
quency behaviours of 3?e 7(0;) and Teff . (uj) , assumed to be respectively proportional 
to \u;\^-^ (0 < 5 < 2) and to |c^|" (a < 0). 

Our main new results follows. In an aging medium, the velocity cor- 

relation function is different from its equilibrium form. The regression theorem is 
not valid. The anomalous diffusion exponent is given hy u = 5 — a. Since a < 0, 
the diffusion exponent is enhanced as compared to the value it would have, with 
the same 3^67(0;), in a thermal bath. This is due to the fact that the weight of 
slow modes in the noise is increased. 

Several questions remain to be solved. In particular, it will be interesting to 
extend the present study so as to obtain a complete analytic description of the 
velocity correlation function and of the mean square displacement at any times, 
i.e. not only in the asymptotic regime. We have previously solved this question 
at equilibrium, in which case we have shown that Mittag-Leffler functions play a 
central role in this description [8] . The out of equilibrium extension is left out for 
future work. 
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